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We study spin motive forces, i.e., spin-dependent forces, and voltages induced by time-dependent 
magnetization textures, for moving magnetic vortices and domain walls. First, we consider the 
voltage generated by a one-dimensional field-driven domain wall. Next, we perform detailed calcu- 
lations on field-driven vortex domain walls. We find that the results for the voltage as a function of 
magnetic field differ between the one-dimensional and vortex domain wall. For the experimentally 
relevant case of a vortex domain wall, the dependence of voltage on field around Walker breakdown 
depends qualitatively on the ratio of the so-called /3- parameter to the Gilbert damping constant, and 
thus provides a way to determine this ratio experimentally. We also consider vortices on a magnetic 
disk in the presence of an AC magnetic field. In this case, the phase difference between field and 
voltage on the edge is determined by the /3 parameter, providing another experimental method to 
determine this quantity. 

PACS numbers: 75.60.Ch, 72.15.Gd, 72.25.Ba 



I. INTRODUCTION 

One of the recent developments in spintronics is the 
study of spin motive forced and spin pumpinj^l. These 
effects lead to the generation of charge and spin currents 
due to time-dependent magnetization textures. The idea 
of spin motive forces due to domain walls is easily un- 
derstood on an intuitive level: if an applied current in- 
duces domain-wall motion,l^l2l Onsager's reciprocity the- 
orem tells us that a moving domain wall will induce a 
current. This idea was already put forward in the eight- 
ies by Berger^. In the case of a domain wall driven by 
large magnetic fields [i.e., well above the so-called Walker 
breakdown field), a fairly simple approach to the prob- 
lem is justified where one goes to a frame of reference in 
which the spin quantization axis follows the magnetiza- 
tion texture.'^ This transformation gives rise to a vec- 
tor potential from which effective electric and magnetic 
fields are derived. Experimentally, the domain-wall in- 
duced voltage has recently been measured above Walker 
breakdowrJi^', and the results are consistent with this ap- 
proach. It has also been shown that the induced voltage 
well above Walker breakdown is determined from a topo- 
logical argument that follows from the properties of the 
above-mentioned vector potential.^ 

The above approach only captures the reactive con- 
tribution to the spin-motive forces. When the velocity 
of the domain wall is below or just above Walker break- 
down, a theory is needed that includes more contribu- 
tions to the spin motive forces. Renewed interest has 
shed light on the non-adiabatic and dissipative contribu- 
tions to the spin motive forcePHUlthat are important in 
this regime. In this paper, we study this regime. 

The article is organized as follows. In Section |ll] we 
summarize earlier results that give a general framework 
to compute electrochemical potentials for given time- 



dependent magnetization textures. In Sect ion [Till we con- 
sider an analytical model for a one-dimensional domain 
wall and numerically determine the form of the spin accu- 
mulation and the electrochemical potential. The results 
agree with the known results for the potential difference 
induced by a moving one-dimensional domain wall.^'* In 
Section|lV]we turn to two-dimensional systems and study 
a vortex domain wall in a permalloy strip. We use a 
micro-magnetic simulator to obtain the magnetization 
dynamics, and numerically evaluate the reactive and dis- 
sipative contributions to the voltage below and just above 
Walker breakdown and compare with experiment.'^ An- 
other example of a two-dimensional system is a vortex 
on a disk which we treat in Section [V] For small enough 
disks, the magnetic configuration is a vortex. Both ex- 
perimentally and theoretically, it has been shown that 
a vortex driven by an oscillating magnetic field will ro- 
tate around its equilibrium position.tI2'22j This gives rise 
to a voltage difference between the disk edge and cen- 
ter as was recently discussed by Ohe et al. [H]. Here, 
we extend this study by including both the reactive and 
the dissipative contributions to the voltage, that turn out 
to have a relative phase difference. This gives rise to a 
phase difference between the drive field and voltage that 
is determined by the so-called /3-parameter. 



II. MODEL 

The spin-motive force field 'P{x) induced by a time- 
dependent magnetization texture that is characterized 
at position x by a uni t-vector magnetization direction 
m(a;, t) is given bjff^lSl 

= ^ [m • {dtm X V,m) + /3(9tm • V,m)] . (1) 
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This force field acts in this form on the majority spins, 
and with opposite sign on minority spins. In this ex- 
pression, the first term is the well-known reactive termfil 
The seco nd ter m describes dissipative effects due to spin 
relaxatioiP^mi and is proportional to the phenomeno- 
logical /3-parameter, which plays an important role in 
current-induced domain-wall motionl^Hs! The spin accu- 
mulation /is in the system follows fronJSl 



A? 



VVs = - V • F 



(2) 



where Asd = \/tD is the spin-diffusion length, with r 
a characteristic spin-flip time and D the effective spin- 
diffusion constant. Here, we assume that the spin- 
relaxation time is much smaller than the timescale for 
magnetization dynamics. The total electrochemical po- 
tential [i that is generated by the spin accumulation due 
to a non-zero current polarization in the system is com- 
puted frorrfl^ 



(3) 



where the current polarization is given by 7^ = (cr^ — 
cr^)/(CT^ + (T^). Note that there is no charge accumulation 
for a-^ — cr^, with a^(cr^) the conductivity of the majority 
(minority) spin electrons. 

The magnetization dynamics is found from 
Landau-Lifschitz-Gilbert (LLG) equation given by 
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Here, iJnmJm] is the micromagnetic energy functional 
that includes exchange interaction, anisotropy, and ex- 
ternal field, and a is the Gilbert damping constant. 



move a domain wall in the direction Qx. Here, we choose 
(5=1. The domain wall width is indicated by A. 

To study the time-evolution of a domain-wall, we let 
(t){x) — 00 (0 so that the wall is described by time- 
dependent collective coordinates 0o(^)}i called 
position and chirality, respectively. For external fields 
smaller than the Walker-breakdown field, there is no 
domain- wall precession {i.e., the chirality is constant), 
and the domain wall velocity v is constant so that 
dtm = —vdx'm.. Since dyXn = 9zm = 0, we immediately 
see that the first term on the right-hand side of Eq. ([l]) 
vanishes, and that the force is pointing along the x-axis. 
We then find that = {Pvh/2)/{\^ cosh[a;/A]2). Due to 
symmetry we have that dy^ = dz^i = dyHs = S^/ig = 0. 
In Figs. [T] and [2] we plot the spin accumulation and the 
electrochemical potential as a function of x. 
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FIG. 1: (color online) Spin accumulation as a function of 
position for several values of the spin-diffusion length. The 
black dotted line gives the value of the source term. The spin 
accumulation tends to zero for x — > ±oo. 
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III. ONE-DIMENSIONAL DOMAIN WALL 

For one-dimensional problems the voltage difference 
can be easily found. For example, an analytic expression 
for the electric current (which is the open-circuit equiva- 
lent of the chemical-potential difference) was obtained by 
one of us for an analytical model for a one-dimensional 
driven domain walLl^ In this section, we solve the po- 
tential problem for a one-dimensional domain wall and 
obtain the explicit position dependence of the spin accu- 
mulation and the chemical potential. 

A one-dimensional domain wall {dytn — (9zm = 0) is 
described 



(x, t) = 2arctan |f 



Q[x-X{t)]/\ I 



(5) 



with m — (sin cos 0, sin 6* sin(/), cos 0). Here, Q — ±1 is 
called the topological charge of the domain wall since it 
indicates the way in which an external field affects the 
domain- wall motion, i.e., a field in the direction -t-z will 




FIG. 2: (color online) Electrochemical potential as a function 
of position. Note that the potential is proportional to the 
polarization and that on the horizontal axis the position x is 
multiplied by the spin-diffusion length. 



From Fig. [2] we see that the total potential differ- 
ence A/t = /i(a; — > oo) — fi{x — — oo) is independent of 
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the spin-diffusion length and linear in the parameter /?: 
A/z = hVfiv/X. Note that this result is only valid below 
Walker Breakdown. 

To find the voltage for all fields B we generalize the 
results for the voltage difference in Ref. 13^ for general 
domain- wall charge Q. A general expression for the volt- 
age in one dimension is given bjEH 

2 e J 



We insert the ansatz [Eq. 
Eq. (|6| and find 



(6) 

([5| with 4>{x) = 4>o{t)] into 



A/i = 



'2H 



(7) 



To find a time-averaged value for the voltage we consider 
the equations of motion for a domain wall that is driven 
by a transverse magnetic field lj221\Mi^ contributing to the 
energy — (/B • m, with g > 0. The equations of motion 
for X{t) and 4>Q{t) are ultimately derived from the LLG 
equation in Eq. Q, and given by 

h 



(1 

(1- 



A 



a^sin(2(/)o) , 
- (3^sin(2(/)o) 



(8) 



Here, K± is the out-of-plane anisotropy constant. These 
equations are solved by 




(9) 



where (..) denotes a time average. It follows that the 
voltage difference for general topological charge is 



A^ = -Sign(B) 



Q nv\pg\B\ 



1 + a2 2|e| 



1 + - I Re 

a 



(10) 



Note that the overall prefactor Sign(i3)(3 makes sense: 
inversion of the magnetic field should have the same re- 
sult as inversion of the topological charge. 

In the above, we used a domain-wall ansatz with mag- 
netization perpendicular to the wire direction. Using 
the topological argument by Yang et al.^ one can show 
that the result is more general and holds also for head- 
to-head and tail-to-tail domain walls. Therefore, for a 
one-dimensional domain wall, the reactive and dissipa- 
tive contributions, i.e., the contributions with and with- 
out /3 in the above expression, to the voltage always have 
opposite sign. 



IV. VORTEX DOMAIN WALL 

For more complicated two-dimensional structures the 
spin-motive force field can have rotation and the simpli- 
fied expression in Eq. (j6]) is no longer valid so that we 
need to treat the full potential problem in Eqs. ([l]|3|. 
Motivated by recent experimental resultJ^ we consider 
in this section the voltage induced by a moving vortex 
domain wall. 

We study the magnetization dynamics using a micro- 
magnetic simulatop^^ from which we obtain the magneti- 
zation m(2;, t) . This simulator solves the LLG equation in 
Eq. \A\i. For comparison with the experiment by Yang et 
al.^^'^we simulate a permalloy sample that has the same 
dimensions as this experiment, i.e. 20nni x 500nni x 
32/im, which is divided in 1 x 128 x 8192 lattice points. 
On this sample, we drive a head-to-head vortex domain 
wall by means of a magnetic field that is pointing from 
right to left, such that the vortex moves from right to left. 
For several field strengths, we obtain the magnetization 
m, and its time-derivative which allows us to compute 
the force field F at each lattice point. Next, we solve 
the matrix problem that is the discrete equivalent to the 
potential problem in Eqs. ([2| and ([3]). For details on this 
calculation, see App. [B] 

We first investigate the velocity of the vortex domain 
wall as a function of the applied field. We use the value 
a = 0.02 for the Gilbert-damping parameter to obtain 
the curve in Fig. [Sj 
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FIG. 3; Velocity of the vortex domain wall as a function of the 
magnetic field strength for a = 0.02. Above Walker break- 
down, the velocity is time-averaged. The line is a guide to the 
eye. 



The decrease in velocity for B = 1.5 mT signals 
Walker breakdown. Indeed, up to fields B = 1.4 niT, the 
vortex moves parallel to the long direction of the sample. 
For B — 1.5 mT, the vortex domain wall motio n is m ore 
complicate and has a perpendicular component .'^^'221 
therefore expect that below Walker Breakdown, just like 
for the one-dimensional domain wall, the vortex domain 
wall only has a dissipative contribution to the voltage. 
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Comparison with the experimental results of Ref. [TT] 
shows that our velocity is roughly a factor 2 higher. 
This might be partly caused by a difference in damping 
and partly by the presence of defects in the experiment 
which causes pinning and therefore a decrease of veloc- 
ity. The exact value of the Walker breakdown field is 
hard to compare, since this depends also on the exact 
value of the anisotropy. Nonetheless our value for the 
Walker breakdown field is of the same order as Ref. [TT] . 
Moreover, what is more important is the dependence of 
wall velocity and wall-induced voltage on the magnetic 
field normalized to the Walker-breakdown field, as these 
results depend less on system details. 
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FIG. 4: Electrochemical potential as a function of position 
for a moving vortex domain wall on the sample. The num- 
bers on the horizontal axes correspond to lattice points with 
separation a = 3.9nm. This specific figure is for a = 0.02, 
H = O.SmT (i.e. below Walker breakdown), P = 1 and 
Asd = a. Note that the peak signals the position of the vortex 
core. 

An example of a specific form of the electrochemical 
potential on the sample due to a field-driven vortex do- 
main wall is depicted in Fig. [4] We see that there is 
a clear voltage drop along the sample, like in the one- 
dimensional model. Additionally, the potential shows 
large gradients around the vortex core and varies along 
the transverse direction of the sample. For each field 
strength, we compute the voltage difference as a func- 
tion of time. For field strengths below Walker Breakdown 
we find that, as expected, only the dissipative term con- 
tributes and the voltage difference rapidly approaches a 
constant value in time. This is understood from the fact 
that in this regime, the wall velocity is constant after a 
short time. The dissipative contribution to the voltage 
is closely related to the velocity along the sample, as can 
be seen in Fig. [5] 

Above Walker breakdown the reactive term con- 
tributes. We find that for /3 = a, the oscillations in 
the reactive component compensate for the oscillations 
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FIG. 5: (color online) reactive (blue squares) and dissipative 
(red triangles) contributions to the voltage as a function of 
time. The numbers on the horizontal axis correspond to time 
steps of 0.565 ns. The green line gives the velocity along the 
sample, it is scaled to show the correlation with the voltage. 
These curves are taken for a — 0.02, /3 = a and field strength 
B = 1.6mT. 

in the dissipative component. If wc look closely to Fig. [5] 
we see that the length of the periods is not exactly equal. 
The periods correspond to a vortex moving to the upper 
edge of the sample, or to the lower edge. The difference 
is due to the initial conditions of our simulation. We 
average the voltage difference over time to arrive at the 
result in Fig. [6j We see that the dissipative contribu- 
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FIG. 6: (color online) Voltage drop along the sample for a = 
0.02 and /3 = a. 

tion becomes smaller for fields larger than the Walker 
breakdown field, whereas the reactive contribution has 
the same sign and increases. In fact, for /3 = a, the re- 
duction of the dissipative contribution is exactly compen- 
sated by the reactive contribution. The /? dependence is 
illustrated in Fig. [7j The behavior is fundamentally dif- 
ferent from the one-dimensional domain-wall situation: 
for the vortex domain wall, the dissipative contribution 
has the same sign as the reactive contribution. 

In order to understand the relative sign, we now dis- 
cuss general vortex domain walls. A single vortex (i.e. 
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FIG. 7: (color online) Voltage drop along the sample for a 
0.02 and several values for /3. 



with vorticity q = +1) is described by two parameters: 
the charge p = ±1 indicates wether the central magnetic 
moment points in the positive or negative z direction 
and the chirality C indicates wether the magnetic mo- 
ments align in a clockwise (C — —1) or anti-clockwise 
(C = -1-1) fashion. We have a vortex that is oriented 
clockwise C = — 1. The relative sign is explained from 
a naive computation of the voltage above Walker break- 
down that does not take into account the rotation of the 
spin-motive force field 

Afi (X — J dx [pdtui ■ dxin + m • {dtm x Sa^m)] 



I3v.. 



cx {^S + l)Vx 



(11) 



where 5 is a positive number and we used that above 
Walker breakdown m ~ m(x — Vxt, y — Vyt) with Vx 
and Vy 0. Note that if Vy = (below Walker 
breakdown) the reactive term, i.e., the second term in 
the above expression, indeed vanishes. We used in the 
last line that above Walker breakdown Vx oc —pvy and 
J dxm • (9j,m x dxm) cx —p. The former equality is 
understood from a geometric consideration: consider a 
sample with a vortex characterized by C = 1, p = 1 and 
VxVy < 0. By symmetry, this is equivalent to C = —1, 
p = — 1 and VxVy > 0. It is therefore clear that the sign 
of VxVy depends on either the polarization, or the hand- 
edness of the vortex. Since we know from the vortex 
domain wall dynamics that reversal of the polarization 
reverses the perpendicular velocity,^ we conclude that 
VxVy does not depend on the handedness of the vortex. 
The latter equality is understood from a similar argu- 
ment: // dxdym ■ (9j,m x 9j,m) changes sign under the 
transformation m — > — m. During this transformation, 
both p ^ —p and C — > — C, and therefore their prod- 
uct cannot account for the total sign reversal. Therefore, 
the integral depends on the polarizatiorP^ but not on 
the handedness of the vortex. The positive number S is 
obtained from our numerical simulation, which suggests 



that the magnetic-field dependence of the voltage is 

A/^ = /3B X constant + (1 - /?/a)|A/Xi.cactivc| ■ (12) 

Note that the sign of the relative contributions can also 
be obtained using the topological argument by Yang et 
al^, which gives the same result. 

We compare our results in Fig. [7] with the experiment 
by Yang et al.'^. If we assume that the voltage below 
Walker breakdown lies roughly on the same line as the 
voltages above Walker breakdown, their results suggest 
a slope of lOnV/Oe. For V ~ 0.8, our results suggest a 
slope of (/3/a)14nV/Oe. Taking into account our higher 
velocity, we find that /3 in the experiment is somewhat 
larger than a. The decrease in slope of the voltage in 
Ref. [11] as Walker breakdown is approached from above 
also suggests (3 > a. 

In conclusion, the behavior of the voltage around 
Walker breakdown allows us to determine the ratio /3/a. 
In experiment, the potential difference as a function of 
the applied magnetic field would show an upturn or 
downturn around Walker breakdown as in Fig. [Tj which 
corresponds to /3 < a and /3 > a, respectively. 



V. MAGNETIC VORTEX ON A DISK 

On small disks (of size /im and smaller) of ferromag- 
netic material the lowest energy configuration is a vor- 
tex. It has been shown that one can let the vortex rotate 
around its equilibrium position by applying an AC mag- 
netic field^® ^". This motion gives rise via Eq. Q to a 
spin motive force on the spins, which induces a voltage on 
the edge of the disk relative to a fixed reference voltage, 
e.g. the disk center. Ohe et al^^ have shown that the 
reactive contribution to the spin motive force field can be 
seen as a dipole that is pointing in the radial direction, 
i.e., the divergence of the force field consists of a posi- 
tive and a negative peak along the radial direction (note 
that the divergence of the force field can be seen as an 
effective charge). Rotation of this dipole gives rise to an 
oscillating voltage on the edge of the sample. Here, we 
consider also the dissipative contribution to the voltage. 

We consider a vortex on a disk with radius R that 
moves around its equilibrium position {i.e., the center 
of the disk) at a distance tq from the center of the disk 
with frequency uj. We use as a boundary condition that 
the magnetization on the edge of the disk is pointing 
perpendicular to the radial direction. In equilibrium, the 
micro-magnetic energy density of the form — Jm- V^m — 
K±rn^ is minimized by 



-y 



mx{x, y) = , " cos 
+ y 

my{x,y) 



!arctan(e-'^V^'+'^'/'') 



m^{x,y) =psin 



2 arctan 



2arctan(e-'^V^'+2''/'' 



(13) 
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where the center of the vortex is chosen at x = y = 0. 
Here k = y/ K±/J is the typical width of the vortex core. 
For permalloy this length scale is of the order ^ lOnm. 
The parameters p and C are defined as before, for defi- 
niteness we choose p=l,C = — 1. To describe clockwise 
circular motion of the vortex around its equilibrium po- 
sition at fixed radius rg we substitute a; — > a; — rg sin(wt) 
and y — 7- y — rocos(a;t). Note that we assume that the 
form of the vortex is not changed by the motion, which 
is a good approximation for rg <C R. 



From the magnetization in Eq. (13), we compute the 
force field using Eq. ([!]). The reactive and dissipative 
contributions to the divergence of the force field are 
shown in Fig. |8] The direction of the dipoles follows 



0.1 



^ 0. 




FIG. 8: The reactive (left) and dissipative (right) contribu- 
tions to the divergence of the force field. White means posi- 
tive values, black is negative values. The reactive contribution 
can be seen as a dipole in the radial direction. The dissipative 
contribution is a dipole perpendicular to the radial direction. 
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FIG. 9: (color online) The reactive (red dashed curve) and 
dissipative (blue dashed curve) contributions to the voltage 
difference between opposite points on the edge of the disk. 
The green line gives the total voltage difference A/xtotai = 
A/^rcac + A^diss, lu this example for /? = 0.4. We used ro = 
lOAsd, R = lOOAad and k = Xsd- For realistic spin-diffusion 
length Asd — 5 nm, these parameter values agree with the 
system of Ohe et al^^. 



structure is lost. 

The dissipative contribution becomes important for 
large values of /3. In principle, it is possible to determine 
/3 by looking at the shift of the peak in the total volt- 
age with respect to the peak in the reactive contribution, 
which is in turn determined by the phase of the applied 
magnetic field. The phase difference between applied 
field and measured voltege then behaves as tan(A(/)) oc /3. 



directly from Eq. ([T]) if we realize that for our system 
—dtm ■ Vm = v{dym)^ is always pointing in the direc- 
tion of the velocity which shows that the dissipative con- 
tribution points along the velocity. Likewise the reactive 
contribution is always pointing perpendicular to the ve- 
locity. 

From the relative orientations of the effective dipoles, 
we expect that the peaks in the reactive and dissipative 
contributions to the voltage on the edge will differ by a 
phase of approximately tt/2 (for rg/i? — this is exact). 
We divide our sample in 1000 rings and 100 angles and 
use the general method in App. |B]to find the voltage on 
the edge shown in Fig. [oj To compare with Ref. pi], 
we take a frequency w/(27r) — 300 MHz and V = 0.8, 
which yields amplitudes for the reactive contribution of 
~ on the edge. However, Ohe et al. suggest that 
voltage probes that are being placed closer to the vortex 
core measure a higher voltage. This indeed increases the 
voltage up to order ~ 10 at r = 2rg. Placing the 
leads much closer to the vortex core does not seem to 
be realistic because of the size of the vortex. Since the 
voltage scales with velocity, it can also be increased by a 
larger radius of rotation, i.e. by applying larger magnetic 
fields. However, for disks larger than 1 /im, the vortex 



VI. DISCUSSION AND CONCLUSION 

We have investigated the voltage that is induced by a 
field-driven vortex domain wall in detail. In contrast to 
a one-dimensional model of a domain wall, the reactive 
and dissipative contribution to the voltage have the same 
sign. The qualitative differences for different values of /3 
provide a way to determine the ratio /3/a experimentally 
by measuring the wall-induced voltage as a function of 
magnetic field. To this end the experimental results in 
Ref. [TT] are in the near future hopefully extended to 
fields below Walker breakdown, which is challenging as 
the voltages become smaller with smaller field. 

We also studied a magnetic vortex on a disk. When 
the vortex undergoes a circular motion, a voltage is in- 
duced in the sample. Earlier work computed the reactive 
voltage on the edge of the disk,'2ll here we include also the 
dissipative contribution to the voltage. We find that the 
phase difference between voltage and AC driving field is 
determined by the /3-parameter. 
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Appendix A: Boundary conditions 



As a boundary condition for the potential problems, 
we demand that the total spin current and charge current 
perpendicular to the upper and lower boundaries is zero: 
= - jJ- = and = + = 0. Therefore, 
the majority-and minority spin currents are necessarily 
zero. They are given by = a^{F — V fj,^) and — 
(Jl{—F — V^^). From this, the boundary conditions on 
the derivatives of the potentials follow as d±iJ,s = d±{fi^— 
/Lt^)/2 = F and d±fj. = d±{iJ,^ + /i^)/2 = 0. We consider 
a two-dimensional sample that is infinitely long in the 
x-direction, and of finite size 2A in the y-direction the 
boundary conditions are 



dyfis{x,y = ±A) = Fy{x,y = ±A) 



(Al) 



To measure the induced voltage, we also put the deriva- 
tives of the potential at infinity to zero so that the bound- 
ary conditions for the electrochemical potential are 

dyii{x,y = ±A) = , 

^ ±c5o, = . (A2) 

Appendix B: Potential problem on a Lattice 

We consider a two-dimensional lattice, where we have 
spin accumulation and an electrochemical potential 
/i*'^ at site Between sites and + 1), there 
can be a particle current density of majority spins 



»J + l/2 



f: 



J + l/2 A^t 



F 



(Bl) 



with aj'"'^^^^ the lattice spacing in the j direction between 
sites and and a particle current density of 

minority spins 



*J + l/2 



a, -F 



7+1/2 A*^ ^^ 



i,J + l/2 



„ f pij + 1/2 c , ij + 1/2 



(B2) 



and equivalently for currents in the i direction. The 
derivative Sj is defined as S^O^'^ — (0*J+i/2 — 
(92j-i/2^^Q«J^ and likewise for (5;. Note that upper in- 
dices (i, j) denote a position on the lattices and lower in- 
dices i or j denote a direction. We can write fi^ = fi + fig 



and /i^ = fi — fj,s- The continuity- like equations for the 
density of majority- and minority spins are (note that 
spins move in the direction of the current) 



(B3) 



with characteristic spin-fiip time r and with the dimen- 
sionless operator A given by 

(B4) 

These definitions allows for non-square lattices with sides 
at position (i ± 1/2, j) or ± 1/2) that have length 
^i±i/2j ^jj±i/2 j-jQ^gj, index denotes the normal di- 
rection), respectively, and the area of the site itself given 
by A''^. 

The equation for the electrochemical potential is ob- 
tained from the continuity equation 

Tl — [- Tl 

(at + a^)A{r'-'"[-V'-'' + V{F''^' - . 



(B5) 



where the current polarization is given by 7^ = {a^ — 
CT^)/(crt + <J^). This result was already obtained for a 
continuous system in Ref. |14) . To find an equation for 
the spin accumulation, we write 

i.j i,j 

a^A[P'^{F'^= - 5^^^)] - a;A[f ^^(-F*'^ - 5^il')] = 

(at + a^)A{t'^[F''^ - S^il'^ - rS^l'•^]} = 

(at + a^)(l - V')A[t^^iF^-'^ - S^i^^)] . (B6) 

If we compare this in the case of a square lattice to the 
expression in ReflHl 



1 

X2~ 



(B7) 



we find that the density of spins that pile up can be 
expressed in terms of the spin accumulation as (n|'-' — 

nl^)/T = (at + a;)(l - V^)fil'3 /{\e\X'^,^). We insert this 
expression to find that the spin accumulation on a lattice 
is determined by 



^ = ±-A[e'^{F^'^-S^^i:^)]. (B8) 
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